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[i]  This  paper  (1)  compares  the  bottom  flows  of  three  existing  high-resolution  global 
simulations  of  the  oceanic  general  circulation  to  near-bottom  flows  in  a  current  meter 
database  and  (2)  estimates,  from  the  simulations,  the  global  energy  dissipation  rate  of  the 
general  circulation  by  quadratic  bottom  boundary  layer  drag.  The  study  utilizes  a 
data-assimilative  run  of  the  Naval  Research  Laboratory  Layered  Ocean  Model  (NLOM),  a 
nonassimilative  run  of  NLOM,  and  a  nonassimilative  run  of  the  Parallel  Ocean 
Program  z-level  ocean  model.  Generally  speaking,  the  simulations  have  some  difficulty 
matching  the  flows  in  individual  current  meter  records.  However,  averages  of  model 
values  of  |u/,|3  (the  time  average  of  the  cube  of  bottom  velocity,  which  is  proportional 
to  the  dissipation  rate)  computed  over  all  the  current  meter  sites  agree  to  within  a  factor  of 
2.7  or  better  with  averages  computed  from  the  current  meters,  at  least  in  certain  depth 
ranges.  The  models  therefore  likely  provide  reasonable  order-of-magnitude  estimates  of 
areally  integrated  dissipation  by  bottom  drag.  Global  dissipation  rates  range  from  0.14 
to  0.65  TW,  suggesting  that  bottom  drag  represents  a  substantial  sink  of  the  ~1  TW 
wind-power  transformed  into  geostrophic  motions. 

Citation:  Arbic,  B.  K.,  J.  F.  Shrivcr,  P.  J.  Hogan,  H.  E.  Hurlburt,  J.  L.  McClean,  E.  J.  Metzger,  R.  B.  Scoti,  A.  Sen,  O.  M.  Smedstad, 
and  A.  J  Wallcraft  (2009),  Estimates  of  bottom  flows  and  bottom  boundary  layer  dissipation  of  the  oceanic  general  circulation  from 
global  high-resolution  models,  J.  Geophys.  Res .,  114 ,  C02024,  doi:10.1029/2008JC005072. 


1.  Introduction 

[2]  The  first  goal  of  this  paper  is  to  eomparc  the  bottom 
flows  in  high-resolution  oeean  general  eireulation  models  to 
flows  in  near-bottom  moored  eurrent  meters.  Beeause  of 
improvements  in  computer  power,  bathymetric  products, 
surface  forcing  produets,  numerical  sehemes,  and  subgnd- 
scale  parameterizations,  the  realism  of  high-resolution 
ocean  models  has  continued  to  improve  in  recent  years. 
Though  some  comparisons  of  high-resolution  models  to 
data  have  focused  on  the  full  water  column  [e.g.,  Penduff  et 
al ,  2006],  most  sueh  comparisons  have  foeused  on  the 
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surfaec  and  upper  oecan  [c.g.,  McClean  et  a l .,  2002; 
Maltrud  and  McClean ,  2005].  Here  we  eonduct  a  compre¬ 
hensive  comparison  of  bottom  flows  in  global  models  with 
those  observed  in  current  meters. 

[3]  The  emphasis  on  bottom  flows  is  motivated  by  the 
second  goal  of  the  paper,  which  is  to  estimate,  from  the 
models,  the  energy  dissipation  rate  of  the  oceanic  general 
circulation  by  bottom  boundary  layer  drag.  In  recent  years 
considerable  effort  has  been  expended  in  quantifying  the 
sources  and  sinks  of  mechanical  energy  for  the  oceanic 
general  eireulation  [e.g.,  Munk  and  Wunsch ,  1998;  Wunsch 
and  Ferrari ,  2004].  One  important  energy  source  is  the 
'*'-'0.85-1  TW  wind-power  input  transformed  into  geo- 
strophie  flows  [e.g.,  Wunsch ,  1998;  Scott,  1999,  Huang  et 
al ,  2006;  Scott  and  Xu,  2008],  How  this  energy  input 
ultimately  dissipates  is  not  well  understood.  Ferrari  and 
Wunsch  [2008]  ean  be  consulted  for  a  detailed  overview. 
Candidate  dissipation  meehanisms  include  internal  wave 
generation  and  breaking  instigated  by  geostrophic  flows 
over  rough  topography  [eg.,  Naveira-Garahato  et  al ., 
2004;  Marshall  and  Naveira-Garahato,  2008;  Nikurashin, 
2008]  and  energy  transfer  from  geostrophic  motions  into 
internal  waves  and  other  submesoscale  motions  in  the  upper 
ocean  [c.g.,  Muller  et  al,  2005;  Polzin,  2008].  Another 
candidate  mechanism  is  bottom  boundary  layer  drag.  We 
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will  focus  on  this  mechanism  here,  though  we  do  not 
discount  the  potential  importance  of  the  other  mechanisms. 

[4]  Several  lines  of  research  point  to  the  importance  of 
bottom  boundary  layer  drag  for  the  dynamics  and  energy 
budget  of  the  oceanic  general  circulation  and  its  associated 
mesoseale  eddy  field.  Weatherly  [1984]  estimated  energy 
dissipation  rates  in  bottom  boundary  layers  for  the  Gulf 
Stream  region,  and  found  that  they  may  balance  nearly  all  of 
the  local  wind-power  input.  A  series  of  papers  [Arbic  and 
Flierl ,  2003,  2004;  Arbic  et  al. ,  2007;  Arbic  and  Scott , 
2008]  have  focused  on  the  sensitivity  of  idealized  geo- 
strophie  turbulence  models  to  bottom  drag,  and  argued  that 
the  horizontal  scales,  vertical  structure,  and  amplitude  of  the 
modeled  eddy  kinetic  energy  compare  well  to  those  in 
observations  only  when  bottom  drag  is  moderately  strong. 
This  is  true  for  both  linear  and  quadratic  bottom  drag  [Arbic 
and  Scott ,  2008].  Other  studies  [e.g.,  Thompson  and 
Schmitz ,  1989,  and  references  therein;  Smith  et  al. ,  2002; 
Riviere  et  a l 2004;  Cessi  et  al. ,  2006;  Thompson  and 
Young ,  2006,  2007]  have  also  found  a  significant  sensitivity 
of  ocean  general  circulation  models  and  geostrophie  turbu¬ 
lence  models  to  the  strength  of  bottom  drag.  Further 
motivation  for  the  importance  of  bottom  boundary  layer 
drag  comes  from  studies  of  the  ocean  tides.  About  2/3  of  the 
globally  integrated  tidal  dissipation  takes  place  in  shallow' 
seas,  where  tidal  flows  are  exceptionally  strong,  primarily 
by  bottom  boundary  layer  drag  [Egbert  and  Ray ,  2003,  and 
references  therein] 

[5]  Sen  et  al.  [2008]  estimated  the  global  energy  dissipa¬ 
tion  rate  of  oceanic  low-frequeney  flows  by  quadratic 
bottom  boundary  layer  drag  (hereafter,  often  referred  to 
simply  as  the  ‘'dissipation  rate”)  using  observations  alone 
The  dissipation  estimates  ranged  from  0.21  to  0.83  TW  (in 
waters  deeper  than  3000  m),  a  substantial  fraction  of  the 
wind-power  input.  The  dissipation  was  concentrated  in 
regions  of  strong  currents  such  as  the  Gulf  Stream,  consis¬ 
tent  with  the  regional  analysis  of  Weatherly  [1984].  Sen  et 
al.  [2008]  utilized  a  database  of  moored  near-bottom  eurrent 
meters.  Relationships  computed  at  the  mooring  sites  be¬ 
tween  surface  flows  measured  by  satellite  altimetry  and 
bottom  flows  measured  by  eurrent  meters  were  exploited  to 
infer  global  maps  of  bottom  flows  from  the  global  satellite 
measurements  of  surface  flows.  Here  we  compute  the 
dissipation  rate  from  direct  estimates  of  bottom  flow,  taken 
from  the  high-resolution  models. 

[6]  We  use  models  with  differing  horizontal  resolutions, 
vertical  resolutions,  and  vertical  coordinate  formulations. 
We  examine  both  data-assimilative  and  nonassimilative 
simulations  of  the  Naval  Research  Laboratory  Layered 
Ocean  Model  (NLOM)  [Shriver  et  al .,  2007],  as  well  as  a 
nonassimilative  simulation  of  the  Parallel  Ocean  Program 
(POP)  [Maltrud  and  McClean ,  2005].  Both  of  the  NLOM 
simulations  are  run  at  1/32°  horizontal  resolution,  while 
POP  is  run  at  1/10°  horizontal  resolution.  The  greater 
horizontal  resolution  of  the  NLOM  simulations  is  potential¬ 
ly  an  important  advantage  for  the  modeling  of  bottom  flows. 
Hurlburt  and  Hogan  [2000]  and  Hogan  and  Hurlburt 
[2000]  showed  that  the  abyssal  eddy  kinetic  energy  in  1/ 
32°  North  Atlantic  and  Japan/East  Sea  regional  simulations 
of  NLOM  was  little  different  from  that  in  1/64°  simulations. 
However,  in  1/16°  runs,  abyssal  eddy  kinetic  energy  was 
notieably  lower,  and  in  1/8°  runs  it  was  much  lower.  Thus  at 


least  in  NLOM,  1/32°  resolution  is  required  to  achieve 
convergence  for  abyssal  eddy  kinetic  energy.  (See  Wallcraft 
et  al.  [2005]  for  a  study  of  the  numerical  convergence  of 
NLOM  with  respect  to  values  of  the  eddy  viseosity.) 
Consistent  with  these  NLOM  results,  Penduff  et  al.  [2006] 
showed  that  the  abyssal  eddy  kinetic  energies  in  1/6° 
simulations  of  the  CLIPPER  model  of  the  Atlantic  Ocean 
were  substantially  weaker  than  those  recorded  by  eurrent 
meters.  On  the  other  hand,  POP  has  greater  vertical  resolu¬ 
tion  than  NLOM  (40  z-levels  versus  6  Lagrangian  layers), 
and  a  more  standard  treatment  of  bottom  topography  (full- 
cell  topography  in  a  z-level  model).  As  described  in  section 
2.3,  the  topography  used  in  NLOM  is  compressed  in  a 
somewhat  nonstandard  way.  Vertical  resolution  and  the 
treatment  of  topography  thus  represent  an  important  poten¬ 
tial  advantage  for  POP. 

[7]  Prior  to  undertaking  our  analysis  on  a  global  scale,  wc 
will  analyze  the  North  Atlantic  sectors  of  the  three  global 
models.  The  density  of  near-bottom  eurrent  meters  we  can 
compare  our  models  to  is  higher  in  the  North  Atlantic  than 
in  other  regions.  Furthermore,  in  a  North  Atlantic  analysis 
we  can  include  a  fourth  high-resolution  simulation,  the 
North  Atlantic  POP  run  of  Smith  et  al.  [2000].  As  noted 
by  Maltmd  and  McClean  [2005],  North  Atlantic  POP 
simulates  the  Gulf  Stream  more  accurately  than  does  global 
POP  As  we  will  see.  North  Atlantic  POP  compares  more 
closely  to  individual  eurrent  meter  records  than  the  other 
three  simulations,  and  thus  serves  as  a  useful  point  of 
comparison. 

[k]  We  reiterate  that  this  paper  utilizes  the  output,  such  as 
it  is,  of  existing  high-resolution  global  simulations.  We  do 
not  undertake  new  simulations  for  this  paper,  as  such 
simulations  are  extremely  expensive  computationally.  The 
ideal  model  output  for  an  analysis  such  as  ours  with  a  focus 
on  low-frequeney  motions  would  consist  of  either  very 
frequently  sampled  (for  instance,  hourly)  output,  from 
which  high-frequency  motions  can  be  removed  with  a 
low-pass  filter,  or  temporal  averages,  i.e.,  daily  averages. 
Hourly  sampling  of  the  models  is  not  possible,  at  least  for 
all  model  variables  and  all  grid  points  in  the  vertical  and 
horizontal  directions,  because  of  the  unfeasibly  large  stor¬ 
age  and  analysis  computers  that  this  would  require  The 
model  output  we  use  here  consists  of  “snapshots”  rather 
than  temporal  averages,  meaning  that  there  will  be  some 
aliasing  of  high-frequency  motions  into  our  analyses  [Jayne 
and  Tokmakian ,  1 997].  Later  we  will  examine  the  impact  on 
our  results  of  the  relatively  infrequent  temporal  sampling  of 
the  NLOM  and  POP  simulations 

[9]  We  will  not  endeavor  to  compute  complete  energy 
budgets  of  these  models,  as  this  would  be  extremely  time- 
consuming  given  the  vast  size  of  the  model  outputs.  Energy 
budgets  of  NLOM  have  been  computed  over  limited 
regions,  and  with  only  two  Lagrangian  layers  in  the  vertical 
direction  In  these  smaller  computations  it  was  found  that 
the  energy  inputs  and  outputs  did  indeed  balance  well,  with 
eddy  viseosity  and  bottom  drag  together  accounting  for  all 
of  the  dissipation  [Hurlburt  and  Thompson ,  1982;  Hurlburt , 
1986;  Schmitz  and  Thompson ,  1993].  Here  we  will  simply 
compute  the  bottom  drag  dissipation  term,  by  inserting  the 
bottom  flows  from  models  into  standard  formulae  for  the 
dissipation  rate.  We  will  eompute  the  wind-power  input  for 
one  of  the  simulations  (nonassimilative  NLOM),  to  test 
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whether  it  is  comparable  to  the  inputs  computed  from 
observations  in  other  studies. 

[10]  While  the  wind-power  input  is  specifically  calculated 
as  an  energy  source  for  geostrophie  motions,  the  bottom 
flows  in  both  current  meter  data  and  high-resolution  models 
are  not  strictly  geostrophie.  Current  meter  data  can  be 
filtered  in  the  time  domain  to  remove  high-frequency 
ageostrophie  motions  such  as  inertial  motions  and  tides. 
What  remains  has  a  strong  overlap  with  geostrophie 
motions.  For  example,  Wunsch  [1997]  found  a  good  visual 
correspondence  between  maps  of  surface  geostrophie 
motions  computed  from  satellite  altimetry  and  the  surface 
extrapolation  of  low-frequency  motions  computed  via  cur¬ 
rent  meters.  Thus  Sen  et  al.  [2008]  argued  that  it  is 
meaningful  to  compare  the  ~1  TW  power  input  trans¬ 
formed  into  geostrophie  flows  with  the  dissipation  of  low- 
frequeney  flows.  Later  we  will  show  that  the  abyssal  flows 
in  nonassimilative  NLOM,  at  least,  are  predominantly  low' 
frequency. 

2.  Ocean  Models 

2.1.  North  Atlantic  POP 

[n]  We  analyze  year  1999  of  the  North  Atlantic  POP 
simulation  of  Smith  et  al.  [2000].  This  simulation  has  40  z- 
levcls  in  the  vertical  direction.  The  model  output  is  saved 
as  snapshots  every  10  days.  The  horizontal  grid  for  the 
output  that  we  analyzed  (downloaded  from  http://www. 
carthsystemgrid.org)  ranges  from  19.9°S  to  72.6°N  in  the 
north -south  direction,  and  from  98°W  to  17.2°E  in  the 
cast  west  direction.  Our  North  Atlantic  analysis,  of  North 
Atlantic  POP  and  of  the  North  Atlantic  sectors  of  the  three 
global  models,  will  be  performed  over  this  North  Atlantic 
POP  domain.  The  longitudinal  bounds  of  this  domain 
exclude  the  Gulf  of  Mexico.  In  North  Atlantic  POP,  the 
Gulf  of  Mexico  is  inserted  into  Africa  as  a  reentrant  feature 
However,  sinee  the  global  models  are  not  run  in  this  way, 
for  simplicity  we  omit  the  Gulf  of  Mexico  from  the  North 
Atlantic  analysis  of  all  four  models.  The  latitudinal  spacing 
of  North  Atlantic  POP  is  1/10°  at  the  equator,  and  Mercator 
elsewhere,  reducing  to  about  0.03°  at  high  latitudes.  The 
longitudinal  spacing  is  1/10°  everywhere.  Every  vertical 
level  of  the  North  Atlantic  POP  output  has  1 280  by  992  grid 
points  (in  the  north-south  and  east-west  directions, 
respectively). 

[12]  As  described  by  Smith  et  al.  [2000],  bottom  topog¬ 
raphy  for  the  North  Atlantic  POP  simulation  was  derived 
from  the  1/12°  ETOP05  database  [National  Oceanographic 
and  Atmospheric  Administration  NOAA ,  1986].  Following 
interpolation  of  this  data  set  to  the  1/10°  Mercator  grid,  the 
depth  at  each  horizontal  grid  point  was  set  equal  to  that  of 
the  nearest  vertical  level  in  the  model.  Minor  modifications 
were  then  made  at  the  Straits  of  Florida,  Faroe  Bank 
Channel,  and  Strait  of  Gibraltar.  The  wind  stress  was 
derived  from  the  6  hourly  ECMWF  TOGA  Global  Surface 
Analysis,  and  then  linearly  interpolated  to  eaeh  model  time 
step  to  avoid  excitation  of  spurious  inertial  oscillations 
[Jayne  and  Tokmakian ,  1997].  (Note  that  this  temporal 
interpolation  was  done  for  the  global  POP  and  NLOM 
simulations  as  well).  Other  details  of  the  surface  forcing, 
c.g.,  surface  heat  and  freshwater  fluxes,  and  temperature 


and  salinity  restoration  at  the  model  horizontal  boundaries, 
arc  described  by  Smith  et  al.  [2000]. 

2.2.  Global  POP 

[13]  We  analyze  year  2003  of  a  global  POP  simulation,  an 
extension  of  the  multiyear  simulation  described  at  length  by 
Maltmd  and  McClean  [2005].  Our  global  POP  analysis  is 
also  based  on  snapshots  put  out  every  10  days.  Global  POP 
employs  a  displaced  pole  grid  in  the  Northern  Hemisphere 
to  allowr  inclusion  of  the  Arctic  Ocean  The  horizontal  grid 
spacing  of  global  POP  is  1/10°  at  the  equator  In  the 
southern  hemisphere,  the  grid  is  Mercator,  while  the  north¬ 
ern  hemisphere  grid  is  smoothly  distorted  to  allow  for  the 
displaced  pole.  The  logical  dimension  of  one  level  of  global 
POP  is  2400  by  3600. 

[14]  Maltmd  and  McClean  [2005]  should  be  consulted 
for  lull  details  of  the  model  bathymetry  and  forcing  fields, 
but  we  repeat  a  small  portion  of  that  discussion  here. 
Nonpolar,  Arctic,  and  Southern  Ocean  bathymetry  was 
derived  from  Smith  and  Sandwell  [1997],  Jakobsson  et  al. 
[2000]  (International  Bathymetric  Chart  of  the  Arctic 
Oeean),  and  Lythe  and  Vaughan  [2001]  (BEDMAP),  re¬ 
spectively.  Again  full  cell  topography  was  used,  and  again 
the  model  depth  field  was  hand-modified  in  regions  of 
important  sills  and  channels  such  as  the  Canadian  and 
Indonesian  archipelagos,  Denmark  and  Faroe-Shetland  over¬ 
flows,  etc.  Surface  fluxes  were  calculated  using  a  combina¬ 
tion  of  the  NCEP/NCAR  rcanalysis  products  [Kalnay  et  al. , 
1 996],  as  well  as  monthly  data  from  various  sources  listed  by 
Maltmd  and  McClean  [2005].  Wind  stress  was  calculated 
offline  using  the  formulation  of  Large  and  Pond  [1982] 
with  elimatologieal  sea  surface  temperature  taken  from 
Shea  et  al.  [1990].  Temperature  and  salinity  distributions 
were  generated  using  a  combination  of  the  Navy's  1/8°  cli¬ 
matological  product  [Fox  et  al .,  2002]  and  the  Steele  et  al. 
[2001]  data  set  for  the  Arctic.  Subsurface  restoring  was  used 
in  limited  areas  around  some  locations  having  overflows 
which  are  difficult  to  simulate  (i.e.,  the  Mediterranean  Sea, 
Red  Sea,  and  Persian  Gulf). 

2.3.  Global  NLOM 

[15]  We  use  1  full  year  of  output  from  a  global  data- 
assimilative  NLOM  simulation  (denoted  by  DANLOM 
here)  (in  current  use  as  a  U.S.  Navy  operational  model 
[Shriver  et  al 2007])  and  1  year  of  output  from  a  global 
nonassimilative  simulation  (denoted  by  NANLOM  here). 
As  described  by  Shriver  et  al.  [2007],  DANLOM  assim¬ 
ilates  sea  surface  height  from  satellite  altimeters  and  sea 
surface  temperature  from  multichannel  satellite  infrared 
radiometers.  The  horizontal  resolution  of  the  NLOM  simu¬ 
lations  is  1/32°  in  latitude  (ranging  from  72°S  to  65°N,  for  a 
total  of  4384  grid  points  in  the  north  south  direction)  and 
45/1024°  in  longitude  (for  a  total  of  8192  grid  points  in  the 
east-west  direction).  Note  that  1/32°  NLOM  has  essentially 
twice  as  many  grid  points  in  each  horizontal  direction  as 
does  1/10°  global  POP.  The  output  we  require  from  DAN¬ 
LOM  and  NANLOM  (bottom  velocities)  was  not  saved 
over  concurrent  years.  Thus  we  analyze  DANLOM  output 
from  calendar  year  2006  and  NANLOM  output  from  2002. 
Snapshots  of  DANLOM  output  were  saved  daily,  while 
NANLOM  snapshots  were  saved  every  3  days. 
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Source 

*  Compressed  (used  in  NLOM) 


Longitude  (degrees) 


Figure  1*  Illustration  of  the  topography  compression  used 
in  NLOM  (sec  text).  Bathymetry  shown  is  along  28°N  in 
the  Atlantic.  The  source  topography  is  taken  from  ET0P05, 
with  minimum  and  maximum  depths  set  at  200  and  6500  m. 
respectively.  The  compressed  topography  is  the  topography 
actually  used  in  NLOM. 


[t6]  NLOM  has  Lagrangian  layers  in  the  vertical  direc¬ 
tion,  and  is  based  on  the  primitive  equation  model  of 
Hurlburt  and  Thompson  [1980],  but  with  greatly  expanded 
capability  [Wallcraft  et  al,  2003].  It  has  6  dynamical  layers 
and  a  bulk  mixed  layer.  The  source  topography  is  a 
modified  version  of  the  1/12°  ETOP05  data  set  \NOAA, 
1986].  The  maximum  allowable  depth  of  NLOM  is  set  to 
6500  m.  The  topography  is  vertically  compressed  to  confine 
it  to  the  lowest  layer  of  NLOM.  This  is  done  by  multiplying 
the  height  of  the  topography  above  the  6500  m  reference 
level  by  a  factor  of  0.69.  Confining  the  topography  to  the 
lowest  layer  avoids  the  problem  of  running  a  Lagrangian 
layer  model  with  vanishingly  thin  layers.  A  solution  to  this 
problem  exists  [e.g.,  Hallberg  and  Rhines ,  1996],  but  has 
not  been  implemented  in  NLOM.  The  lack  of  vanishingly 
thin  layers  in  NLOM  makes  for  greater  computational 
efficiency,  allowing  it  to  be  run  at  horizontal  resolutions 
not  yet  feasible  in  other  ocean  models.  NLOM  has  lateral 
boundaries  that  follow  the  200  m  isobath  with  a  few 
exceptions,  such  as  the  shallow  straits  around  the  Japan/ 
East  Sea.  Additional  details  arc  provided  by  Shiver  et  al. 
[2007,  and  references  therein].  Figure  1  demonstrates  the 
topography  compression  used  in  NLOM.  Differences  be¬ 
tween  the  topography  actually  used  in  NLOM  and  the 
source  topography  arc  largest  in  shallow  waters.  The 
minimum  depth  of  200  m  in  the  source  topography  yields 
a  minimum  depth  in  NLOM  of  2 1 53  m.  As  we  will  sec  later, 
the  large  differences  between  the  source  depths  and  the 
depths  used  in  NLOM  in  shallow  regions  yield  flows  in 
these  regions  that  arc  weaker  than  those  in  observations. 

[n]  The  NLOM  simulations  were  forced  with  a  hybrid  of 
the  Fleet  Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  Navy  Operational  Global  Atmospheric  Predic¬ 
tion  System  (NOGAPS)  [Rosmond  et  al .,  2002]  and  the 
Hellerman  and  Rosens tein  [1983]  wind  stresses.  The  long¬ 
term  mean  (August  1990  to  July  1999  for  FNMOC)  is 
subtracted  from  the  FNMOC  stresses  and  replaced  by  the 


annual  mean  from  the  Hellerman  and  Rosenstein  [1983] 
data  set. 

3.  Current  Meter  Data 

[18]  As  in  the  paper  by  Sen  et  al.  [2008],  we  use  moored 
current  meter  data  obtained  from  the  Deep  Water  Archive 
and  Buoy  Group  Archive  of  Oregon  State  University’s 
Buoy  Group  (http://cmdac.oce.orst.cdaeds.html,  or  http:// 
emrccords.net).  Since  bottom  boundary  layer  depths  arc  of 
order  10  m  [Weatherly  and  Martin ,  1978],  we  require  the 
current  meters  to  lie  at  least  10  m  above  the  sea  bottom  We 
define  a  current  meter  as  sufficiently  near  the  seafloor  if  it  is 
less  than  10%  of  the  seafloor  depth  above  the  seafloor. 
When  multiple  instruments  from  the  same  mooring  fit  our 
criteria,  only  the  deepest  instrument  is  used  Tides  and  other 
high-frequency  motions  arc  removed  with  a  72-hour  low- 
pass  filter.  In  records  that  contain  data  gaps,  the  filter  is 
applied  separately  to  each  segment  of  complete  data.  Only 
cuiTcnt  meter  records  exceeding  180  days  in  length  are 
retained.  Typical  lengths  of  the  current  meter  records  we  use 
range  from  0.5-2  years.  While  the  current  meter  database 
spans  many  years  (1970s  to  present),  the  model  output  we 
analyze  is  confined  to  specific  years,  listed  in  section  2. 

[19]  We  select  only  those  mooring  locations  for  which  a 
modcl/data  comparison  involving  all  three  global  simula¬ 
tions  can  be  made.  The  number  of  current  meters  in  the 


Depths  less  than  1 000  m 
Depths  between  1 000  and  3000  rr 
Depths  exceeding  3000  m 


(a)  North  Atlantic  sector  current  meters 


300  320  340  360 

(b)  Global  current  meters 


Figure  2.  (a)  Locations  of  187  North  Atlantic  current 
meters  deemed  sufficiently  near  the  seafloor  (sec  text)  and  for 
which  global  DANLOM,  NANLOM,  and  POP  output,  as 
well  as  output  from  the  North  Atlantic  POP  simulation,  are 
available  for  comparison,  (b)  Locations  of  382  moored 
current  meters  used  in  the  global  analysis  of  the  global 
DANLOM,  NANLOM,  and  POP  simulations.  Approximate 
seafloor  depths  of  the  moorings  arc  also  indicated  in  Figure  2. 
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database  which  satisfy  all  of  our  criteria  is  382.  A  subset  of 
these,  numbering  1 87,  can  additionally  be  compared  to  the 
North  Atlantic  POP  simulation.  Figure  2a  shows  the  loca¬ 
tions  of  the  1 87  current  meters  used  in  the  North  Atlantic 
analysis,  while  Figure  2b  shows  the  382  current  meters  used 
in  the  global  analysis. 

4.  Dissipation  Formula 

[20]  The  time-averaged  dissipation  rate  at  a  model  gnd 
point  (or  mooring  location)  is  computed  as 

D(M)  =  PQkP  (1) 

[ Taylor ,  1919],  where  0  and  4>  arc  respectively  the  longitude 
and  latitude  of  the  grid  point,  p  =  1035  kg  m  3  is  the 
average  density  of  seawater,  |u^|  is  the  magnitude  of  the 
bottom  velocity  vector,  and  overbars  denote  time  averaging. 
In  DANLOM  and  NANLOM  we  take  the  bottom  velocity  to 
be  the  velocity  in  the  bottom  (sixth)  layer.  In  North  Atlantic 
and  global  POP  the  bottom  velocity  is  taken  from  the  lowest 
active  level  of  the  40  z-levcls.  This  lowest  level  varies  as  a 
function  of  geographical  location.  Only  “total”  velocities 
are  utilized  in  our  study.  A  time  mean  is  not  removed  from 
either  model  or  mooring  velocities.  In  ocean  models,  the 
quadratic  drag  coefficient  cd  typically  lies  between  0.001 
and  0.003.  Observations  indicate  a  somewhat  wider  range  of 
plausible  values  [c.g.,  Trowbridge  and  Lentz,  1998;  Duncan 
et  ai ,  2003].  The  parameterization  (1)  docs  not  cover  the 
complex  internal  wave  dynamics  arising  in  topographically 
rough  areas,  which  require  an  additional  parameterization 
[c.g.,  Jayne  and  St.  Laurent ,  2001;  Gamer ,  2005]. 

[21]  The  native  values  of  cd  used  in  the  model  simu¬ 
lations  arc  0.003/0.002/0.001225/0.001225  for  DANLOM/ 
NANLOM/global  POP/North  Atlantic  POP.  We  will  present 
dissipation  rates  computed  using  the  native  cd  values,  which 
would  be  the  rates  appropriate  for  an  examination  of  model 
energy  budgets.  Another  approach  we  follow  is  to  regard  the 
models  as  providing  estimates  of  bottom  velocity,  w  hich  can 
be  combined  with  any  reasonable  value  of  cd  to  yield 
dissipation  rates,  irrespective  of  considerations  of  model 
energy  budgets.  We  therefore  also  present  dissipation  rates 
computed  using  a  common  value  of  cd  =  0.0025  for  all  the 
simulations.  The  common  value  makes  for  easier  compari¬ 
son  of  results  from  the  different  simulations,  as  well  as  easier 
comparison  with  the  obscrvationally  based  dissipation  esti¬ 
mates  of  Sen  et  al.  [2008],  in  which  cd  was  also  set  to  0.0025. 

5.  Results 

[22]  The  first  goal  of  this  section  is  to  determine  how 
realistic  the  bottom  flows  in  the  models  arc,  through 
comparison  to  flows  measured  by  near-bottom  moored 
current  meters.  The  second  goal  is  to  discuss  the  dissipation 
rates  computed  from  the  model  bottom  flows.  Values  of 
areally  averaged  dissipation  will  be  given,  and  maps  of  the 
dissipation  rates  will  be  shown.  The  analysis  will  be  done 
first  in  the  North  Atlantic,  and  then  on  a  global  scale. 

[23]  The  model/data  comparison  will  focus  on  values  of 

|ui,|\  since  this  quantity  multiplied  by  p  cd  yields  the 
dissipation  rate.  The  model  and  mooring  |u/,|3  values  will 


be  plotted  alongside  each  other  as  a  function  of  seafloor 
depth.  As  noted  earlier,  there  is  reason  to  believe  that 
NLOM  will  not  perform  as  well  in  shallow  waters  as  in 
abyssal  waters.  Additionally,  scattcrplots  will  be  made  of 

model  versus  moonng  values  of  |u/>|3.  We  display  two 
different  measures  of  the  model  values  at  the  mooring  sites: 
(1)  the  value  of  ju/,|3  interpolated  to  the  mooring  location 
and  (2)  the  model  values,  in  1 0  by  1 0  boxes  centered  on  the 
mooring  locations,  which  come  closest  to  matching  the 
mooring  values.  The  two  measures  arc  denoted  by  “Inter¬ 
polated”  and  “Best”,  respectively,  in  this  paper.  A  third 
measure,  the  mean  of  model  values  over  the  1 0  by  1 0  boxes, 
w'as  also  computed.  The  mean  values  arc  in  most  eases 
similar  to  the  interpolated  values.  Hence  for  the  sake  of 
brevity  we  do  not  further  discuss  the  mean  values  in  this 
paper.  The  best  measure  is  more  forgiving  than  the  inter¬ 
polated  measure,  as  it  compares  well  with  the  data  in  eases 

where  model  |ub|3  fields  contain  values  comparable  to  those 
in  the  moonngs,  but  in  slightly  offset  locations.  In  addition 
to  displaying  model/data  comparisons  at  individual  mooring 
sites,  we  will  compare  averages  of  the  two  model  measures 

of  |u/,|3,  taken  over  many  moonng  locations,  to  averages  of 
the  mooring  values  themselves.  This  exercise  aids  in  gaug¬ 
ing  the  skill  of  the  models  at  predicting  integrals  of  |u/>|3 
(hence  dissipation)  over  large  areas.  In  another  test  of  model 
skill,  we  will  compute  the  best  values  from  model  locations 
sampled  randomly,  to  ascertain  whether  such  sampling 
yields  different  results  than  sampling  the  models  at  the 
actual  moonng  locations. 

5.1.  North  Atlantic  Results 

[24]  We  begin  with  plots  of  | u * | 3  versus  seafloor  depth. 

Figure  3a  shows  the  interpolated  values  of  |u/>|3  in  the  North 
Atlantic  POP  simulation,  and  in  the  North  Atlantic  sectors  of 
the  global  DANLOM,  NANLOM,  and  POP  simulations,  at 
the  187  North  Atlantic  mooring  sites,  as  a  function  of  the 

seafloor  depth  at  the  moonngs.  Mooring  values  of  |u/>|3  arc 
also  plotted.  At  the  two  current  meter  sites  having  seafloor 
depths  less  than  1000  m  (both  located  at  the  Strait  of 

Gibraltar),  values  of  |u/>|3  in  North  Atlantic  POP  arc  some¬ 
what  lower  than  those  in  the  moorings.  Interpolated  |u/>|3 
values  in  the  three  global  models  arc  much  too  weak  at  these 
two  locations.  Note  that  in  this  paper  we  denote  locations 
having  seafloor  depths  less  than  1000  m  as  “Shallow.”  In 
“Intermediate”  depths  (depths  between  1000  and  3000  m), 

none  of  the  three  global  models  have  interpolated  | ] 3 
values  lying  near  the  largest  moonng  values,  while  North 
Atlantic  POP  does  have  values  lying  close  to  this  high  level. 
On  the  other  hand,  for  “Abyssal”  waters  (seafloor  depths 

exceeding  3000  m),  the  interpolated  | u ^  | 3  values  in  all  four 
models  cluster  at  a  level  near  that  seen  in  the  moorings. 

[25]  As  is  to  be  expected,  there  is  a  visually  tighter  fit 
between  models  and  data  in  Figure  3b,  which  displays  the 
best  model  values.  Here  North  Atlantic  and  global  POP 
perform  well  over  all  depths,  with  many  of  the  |u/,|3  values 
lying  nearly  on  top  of  those  seen  in  the  moorings.  The 
DANLOM  and  NANLOM  best  values  continue  to  be  much 
too  small  in  shallow  and  intermediate  depths,  but  lie  at 
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(a)  Interpolated  model  values 
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quantify  this  by  comparing  averages  of  model  and  eurrent 
meter  |u*|3  values  over  many  mooring  sites.  We  define 


7  —  ' 


£Xl  \Ubmodel\j 
i — !  I  ^bmooring  j  | 


(2) 


where  i  is  an  index  of  the  mooring  locations,  M  is  the  total 
number  of  moorings  over  which  7  is  computed,  | u  1  ^  is 
the  measure  of  the  model  |u*|3  values  at  mooring  location  /, 
and  \uhmoonng/  is  the  mooring  value  of  |u*|3  at  loeation  /. 
Table  1  gives  values  of  7  in  the  North  Atlantic  for  the  four 
models  discussed,  for  the  two  measures  of  model  |u/>|3 


(b)  Best  model  values 
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Figure  3.  Values  of  | u * | 3  computed  from  187  moored 
near-bottom  North  Atlantic  eurrent  meters,  alongside  values 
at  the  mooring  locations  from  North  Atlantic  POP  and  the 
North  Atlantic  seetors  of  DANLOM,  NANLOM,  and  global 
POP,  plotted  versus  the  seafloor  depths  of  the  moorings. 

(a)  Model  values  interpolated  to  the  mooring  coordinates. 

(b)  Model  values,  taken  from  1°  by  1°  boxes  centered  on 
the  moorings,  whieh  compare  most  closely  to  the  mooring 
values.  Units  of  |u/>|3  are  s  \ 


approximately  the  same  level  as  the  mooring  values  for 
seafloor  depths  exceeding  3000  m. 

[26]  Next  we  display  seattcrplots  of  model  versus  moor¬ 
ing  values  of  |u/>|3.  In  Figures  4a  and  4b  we  respectively 
show  the  interpolated  and  best  model  | u/> | 3  values  at  the 
mooring  sites,  plotted  against  the  mooring  |u/>f  values.  Extra 
lines  are  drawn  at  a  45°  angle  (all  points  would  lie  along  these 
lines  if  the  models  matehed  the  data  perfectly).  In  Figure  4a, 
the  scatter  is  considerable,  even  with  North  Atlantic  POP  The 
best  values  give  a  tighter  fit,  especially  in  the  ease  of  the  two 
POP  simulations.  In  all  four  simulations,  and  for  both 
interpolated  and  best  measures,  model  values  display  a  weak 
bias  with  respeet  to  the  mooring  values,  as  seen  in  the  greater 
number  of  points  below  the  45°  line  than  above 

[27]  The  indications  thus  far  arc  that  the  models  have 
some  difficulty  matching  individual  eurrent  meter  records. 
However,  we  have  also  seen  visual  indications  that  the 
model  | U/>|3  values  cluster  at  approximately  the  same  level 
as  the  moonng  values,  at  least  in  certain  depth  ranges.  We 


Figure  4.  As  in  Figure  3  exeept  that  model  \uh\}  values 
are  plotted  versus  mooring  |u*f  values. 
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Table  1.  Values  of  7  Computed  Over  187  Moored  Near-Bottom 
Current  Meters  in  the  North  Atlantic3 


1 


Model 

Measure 

Shallow 

Intermediate 

Abyssal 

DANLOM 

Interpolated 

0.000078 

0.025 

0.89 

DANLOM 

Best 

0.00013 

0.051 

0.72 

NANLOM 

Interpolated 

0  000019 

0.026 

0.50 

NANLOM 

Best 

0.000047 

0.055 

0.48 

Global  POP 

Interpolated 

0  069 

0.16 

0.85 

Global  POP 

Best 

0.88 

0.90 

0.77 

North  Allantic 

Inierpolated 

0.37 

0.86 

0.73 

POP 

North  Atlantic 

Best 

0.78 

1.0 

0.78 

POP 

North  Allantic 

Best  random 

0.057 

0.15 

0.32 

POP 

aThe  meanings  of  inlerpolaled,  besl,  and  besl  random  as  measures  of  the 
model  are  explained  in  the  text.  Shallow  denotes  locations  with  moonng 
seafloor  depths  less  than  1000  m  (M  =  2),  intermediate  denotes  depths 
between  1000  and  3000  m  (Xf  48),  and  abyssal  denotes  deplhs  exceeding 
3000  m  (M-  137).  See  equation  (2)  and  surrounding  text  for  definition  of 
7.  M  denotes  number  of  current  meters. 


values,  and  for  three  ranges  (shallow,  intermediate,  abyssal) 
of  seafloor  depths  at  the  mooring  sites.  In  shallow  and 
intermediate  waters,  consistent  with  the  discussion  of 
Figure  3,  the  7  values  for  DANLOM  and  NANLOM  are 
much  less  than  unity,  for  North  Atlantic  POP  are  of  order 
unity  (lying  between  0.37  and  1.0)  for  both  measures,  and 
for  global  POP  are  of  order  unity  for  the  best  measure  only. 
In  abyssal  waters,  the  7  values  for  all  models  and  both  model 
measures  are  of  order  one,  lying  between  0.48  and  0.89. 

[28]  Table  1  contains  a  somewhat  counterintuitive  result 
(some  of  the  best  7  values  arc  lower  than  the  corresponding 
interpolated  7  values).  We  use  a  simple  example  to  illustrate 

how  an  interpolated  |u/>|3  value  could  actually  match  the 
\ub\  value  at  a  single  moonng  more  closely  than  docs  the 

best  |u/>|3  value.  Suppose  that  the  mooring  value  is  ‘*4” 
(arbitrary  units,  for  simplicity)  and  that  the  four  surrounding 
model  grid  points,  all  lying  equal  distances  from  the  moor¬ 
ing,  have  |u/,|3  values  of  1,  3,  6,  and  6.  The  model 
interpolated  |u/>|3  value  would  then  be  4,  exactly  matching 
the  mooring  value,  while  the  model  best  |u/>|3  value  of  3 
does  not.  Because  it  is  possible  at  a  single  mooring  for 
the  interpolated  |u/,|3  value  to  perform  better  than  the  best 

|u 6p  value,  it  is  also  possible  for  interpolated  values  to 
perform  better  when  averages  are  taken  over  many  mooring 
sites. 

[29]  It  is  fair  to  ask  whether  order  one  values  of  7  reall> 
do  demonstrate  model  skill  in  matching  averages  of  |u/,  | 3 
over  large  areas.  This  is  particularly  true  given  that  the  best 
values  are  constructed  to  compare  well  to  the  mooring 
values.  On  the  other  hand,  the  best  values  will  only  compare 
well  to  observations  if  small  subdomains  of  the  model 
contain  at  least  one  grid  point  having  flows  of  approximately 
correct  magnitude.  With  this  in  mind,  we  performed  the 
following  exercise  as  a  further  test  of  model  skill.  For  each 
current  meter  in  a  particular  depth  range,  we  picked  a  grid 
point  lying  in  the  same  depth  range  (as  defined  by  the 
bottom  topography  used  in  the  model),  at  random,  from  the 


North  Atlantic  POP  simulation.  We  then  found  the  model 
|u/>|3  values,  taken  from  1°  by  1°  boxes  centered  on  these 
random  model  grid  points,  which  compare  most  closely  to 

the  mooring  values  of  |u/>|'\  The  exercise  was  repeated 
several  times,  and  an  average  was  taken  over  all  the 
exercises,  to  settle  on  a  value  of  7  denoted  “Best  Random” 
in  the  last  line  of  Table  1 .  The  best  random  7  values  for 
North  Atlantic  POP  are  substantially  smaller  than  the  best 
and  interpolated  values,  for  all  three  depth  ranges.  The  best 
random  values  were  also  calculated  in  the  abyss  for  the  three 
global  models.  The  topography  used  to  define  seafloor 
depths  in  DANLOM  and  NANLOM  for  this  exercise  is 
the  source  topography,  i  e.,  the  topography  before  the 
vertical  compression  used  in  NLOM  takes  place.  The 
abyssal  best  random  values  are  0.77  for  DANLOM,  0.25 
for  NANLOM,  and  0  29  for  global  POP.  Except  in  the  ease 
of  DANLOM,  these  values  are  substantially  smaller  than  the 
best  and  interpolated  7  values,  suggesting  some  degree  of 
model  skill  in  predicting  abyssal  flows.  We  will  return  to 
this  point  in  section  5.2. 

[30]  Since  many  of  the  7  values  in  Table  1  are  of  order 
unity,  at  least  in  specific  depth  ranges,  we  infer  that  all  of  the 
models  may  yield  reasonable  order-of-magnitude  estimates 
of  areally  integrated  dissipation.  Table  2  lists  the  dissipations 
[D]  integrated  over  the  North  Atlantic,  computed  as 

[5]  =  J  ~D(0,  </>)  dA  =  J  pcd~\^dA.  (3) 

where  the  [  ]  operator  represents  an  areal  integration  over 
the  domain  defined  by  the  North  Atlantic  POP  simulation, 
and  dA  is  an  element  of  area.  Values  arc  given  for  shallow, 
intermediate,  and  abyssal  depths.  Again,  the  source 
topography  is  used  to  define  depths  in  DANLOM  and 
NANLOM  for  these  computations.  Total  denotes  an 
integration  over  all  depths.  For  each  simulation,  dissipation 
rates  are  given  for  two  values  of  c(f.  the  native  value  used  in 
the  simulation,  and  a  common  value  of  0.0025.  We  note  that 
the  native  dissipations  arc  within  a  factor  of  2  or  less  of  the 
common  dissipations,  and  concentrate  our  discussion  on  the 
latter,  since  intereompanson  between  different  simulations 
is  easier  in  this  ease.  The  dissipations  in  DANLOM  and 
NANLOM  take  place  primarily  in  the  abyss.  Conversely,  in 


Table  2.  Values  of  Dissipation  Rate  [£>]  in  the  North  Atlantic3 


[/>] 

Model 

 Si.   

Shallow 

lnlermediate 

Abyssal 

Tolal 

DANLOM 

Native 

0.00059 

0.0079 

0.13 

0.14 

NANLOM 

Native 

0.00042 

0.0029 

0.023 

0.026 

Global  POP 

Native 

0.028 

0.010 

0,0095 

0.048 

North  Atlantic 
POP 

Native 

0034 

0.011 

0.016 

0  061 

DANLOM 

Common 

0.00049 

0.0066 

0.11 

0.11 

NANLOM 

Common 

0  00053 

0.0036 

0.029 

0.033 

Global  POP 

Common 

0.058 

0.021 

0.019 

0.099 

North  Atlanlic 
POP 

Common 

0.069 

0.022 

0.033 

0.12 

aTotal  denotes  an  integration  over  all  deplhs.  Computations  are  done  with 
both  native  values  of  cd  (0.003  for  DANLOM,  0.002  for  NANLOM, 
0.001225  for  global  POP,  and  0.001225  for  North  Atlanlic  POP)  and  a 
common  value  of  cd  -  0.0025.  In  TW. 
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Figure  5.  Maps  of  dissipation  rate  (mW  m  2)  in  the  North  Atlantic  sectors  of  the  three  global  models 
(a)  DANLOM,  (b)  NANLOM,  and  (c)  global  POP,  and  (d)  in  North  Atlantic  POP.  Figures  5a,  5b,  5c,  and 
5d  extend  to  72°N.  A  common  value  of  cd  -  0.0025  was  used  to  make  Figures  5a,  5b,  5c,  and  5d. 


both  North  Atlantic  and  global  POP  a  substantial  fraction 
(about  70  80%)  of  the  total  dissipation  takes  place  in 
shallow  and  intermediate  depths.  While  it  is  not  difficult  to 
find  inconsistencies  between  the  simulations,  there  is  also 
some  room  for  optimism.  The  [D]  values  for  global  POP  arc 
within  a  factor  of  1.8  or  less  of  those  in  North  Atlantic  POP, 
over  all  three  depth  ranges.  The  NANLOM  abyssal  [D] 
value  is  very  close  to  that  in  North  Atlantic  POP  and  is 
within  a  factor  of  1.5  of  that  in  global  POP  The  total  [D] 
values  for  the  four  models  lie  within  a  factor  of  about  4  of 
each  other.  To  us  the  discussion  of  Tables  1  and  2  suggests 
that  the  abyssal  and  total  dissipations  in  all  four  simulations, 
and  the  shallow  and  intermediate  dissipations  in  the  two 
POP  simulations,  arc  accurate  to  factors  of  about  4  or  better. 

[31]  We  now  examine  maps  of  the  dissipation  rates.  For 
ease  of  comparison,  Figure  5  displays  the  results  of  all  four 
simulations  mapped  onto  the  same  grid,  which  has  a 
northern  boundary  at  72°N.  Also  for  ease  of  comparison, 
Figures  5a,  5b,  5c,  and  5d  arc  made  using  the  common  cd 
value  of  0.0025.  The  maps  feature  important  similarities,  for 
instance  a  strong  Gulf  Stream  system,  as  well  as  clear 
differences.  Both  POP  simulations  feature  large  dissipation 
rates  along  the  northeast  coast  of  South  America,  in  waters 
between  Great  Britain  and  Iceland,  along  the  cast  coast  of 
Greenland,  and  in  the  Labrador  Sea,  that  arc  not  present  in 


DANLOM  and  NANLOM.  In  the  Gulf  Stream  region,  large 
dissipations  take  place  over  a  greater  areal  extent  in  DAN¬ 
LOM  and  North  Atlantic  POP  than  in  NANLOM  and  global 
POP.  In  DANLOM  there  are  patches  of  large  dissipation 
just  off  of  the  equator  (starting  at  about  ±5°)  which  arc  not 
seen  in  the  other  models.  Some  of  the  tropical  Atlantic 

current  meters  in  our  database  exhibit  |u^|?  values  compa¬ 
rable  to  those  seen  at  ±5°  off  of  the  equator  in  DANLOM, 
but  these  current  meters  lie  within  1  -2°  of  the  equator.  Thus 
it  is  difficult  to  determine  with  the  data  at  hand  whether 
these  patches  arc  real  or  artificial.  The  patches  may  result 
from  the  scheme  used  in  DANLOM  to  propagate  assimi¬ 
lated  information  at  the  surface  to  deep  layers  of  the  model 
[Hurlburt  et  al .,  1990].  The  statistical  inference  is  per¬ 
formed  globally,  but  the  geostrophic  correction  to  velocities 
is  tapered  off  between  5  and  2  degrees.  The  lack  of  a 
geostrophic  adjustment  to  the  currents  along  the  equator 
after  assimilation  of  sea  surface  height  evidently  results  in 
low  values  of  along-cquatonal  dissipation  in  DANLOM  that 
resemble  those  in  the  other  models.  The  patches  of  large 
dissipation  in  near-equatorial  DANLOM  are  in  part  respon¬ 
sible  for  the  larger  [/)]  and  best  random  7  values  seen  in 
abyssal  DANLOM  than  in  the  abyssal  regions  of  the  other 
models 
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(a)  East  Greenland/Labrador 


(b)  Eastern  North  Atlantic 
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(c)  Gulf  Stream 
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Figure  6.  Current  meters  with  large  time-averaged 
currents  (see  text)  in  four  subregions  of  the  North^  Atlantic 
domain.  Values  of  M  (number  of  moorings)  and  Daverage  = 

0.0025  p  |u/,| '  averaged  over  all  of  the  moonngs  (units  in 
mW  m~2),  are  listed  below,  (a)  East  Greenland/Labrador 
Sea_(Af-  35,  Daverage  -  40).  (b)  Eastern  North^Atlantie  (M  = 
27,  Daverage  =  8.9).  (e)  Gulf  Stream  (Af  =  3 1 ,  Daverage  =  6.7). 
(d)  South  America  (M  -  7,  Daverage  -  6.7).  As  in  Figure  2, 
black  crosses  denote  locations  with  seafloor  depths  exceeding 
3000  m,  and  magenta  eireles  denote  locations  with  seafloor 
depths  between  1000  and  3000  m. 


[32]  The  differences  in  the  four  maps  shown  in  Figure  5 
motivate  us  to  examine  how  well  strong  currents  in  specific 
subregions  are  eaptured  in  the  models.  Figure  6  shows  four 
such  subregions  of  the  North  Atlantic  domain.  In  this 
analysis  we  use  only  current  meters  for  which  0.0025  p 
|u/,|3  exceeds  1  mW  m~  2.  Such  current  meters  account  for 

the  bulk  of  the  |u/>|  "  averaged  aeross  all  of  the  eurrent 
meters  in  these  subregions.  We  exelude  the  two  eurrent 
meters  in  the  Strait  of  Gibraltar,  whieh  would  otherwise 
dominate  the  statistics  of  the  Eastern  North  Atlantic 
(Figure  6b),  and  whieh  have  already  been  extensively 
compared  to  the  models,  in  the  “Shallow”  column  of 
Table  1 .  Table  3  displays  the  7  values  computed  in  the  four 
subregions.  For  the  sake  of  brevity,  only  best  7  values  are 
shown  (interpolated  7  values  arc  generally  lower,  at  times 
much  lower).  All  four  simulations  eapture  the  Gulf  Stream 
well,  with  7  values  ranging  from  0.68  to  1.1  All  four 


simulations  perform  reasonably  well  in  the  region  we 
denote  by  “South  America”,  with  7  values  ranging  from 
0.30  to  0.39.  The  strong  currents  along  the  East  Greenland 
eoast  and  Labrador  Sea  (mainly  in  waters  of  intermediate 
depth)  are  eaptured  well  by  the  two  POP  simulations,  but 
not  by  the  two  NLOM  simulations,  as  expected  given  the 
differences  seen  in  Figure  5.  Similarly,  the  strong  currents  in 
the  Eastern  North  Atlantic  (also  mainly  in  waters  of  inter¬ 
mediate  depth)  are  much  better  eaptured  by  POP  than  by 
NLOM.  The  agreement  between  the  POP  simulations  and 
eurrent  meter  reeords  in  these  intermediate-depth  locations 
in  the  North  Atlantic  gives  us  eonfidenee  that  the  strong 
currents  along  continental  edges  may  be  reasonably  aeeu- 
rate  in  POP.  This  will  be  important  in  seetion  5.2,  where  we 
will  see  that  global  POP  has  strong  currents  along  continental 
edges  throughout  the  globe,  usually  in  locations  where  no 
eurrent  meter  data  is  available  for  model  validation. 

5.2.  Global  Results 

[33]  As  in  the  North  Atlantic  analysis,  we  begin  with 
plots  of  model  |u/>|3  values  versus  seafloor  depth  and  versus 
mooring  |u*|3  values.  Figure  7  shows  the  interpolated  and 
best  values  of  |u/> | 3  in  the  three  global  simulations,  at  the 

382  global  mooring  sites,  together  with  the  moonng  |u/>|3 
values,  as  a  function  of  the  seafloor  depth  at  the  moorings 
Figure  7b,  displaying  best  values,  is  visually  tighter  than 
Figure  7a,  displaying  interpolated  values.  In  shallow  and 

intermediate  depths,  the  interpolated  |u/>|3  values  in  all  three 
global  simulations  are  systematically  smaller  than  the  cur¬ 
rent  meter  |u/>| 3  values.  This  is  especially  true  for  the  two 
NLOM  simulations.  The  best  values  of  global  POP  compare 
more  closely  to  observations  in  shallow7  and  intermediate 
depth  waters  than  do  either  the  interpolated  POP  values,  or 
the  best  values  from  DANLOM  and  NANLOM.  In  abyssal 
waters,  for  both  interpolated  and  best  measures,  all  three 

global  simulations  have  |u *>|3  values  that  appear  to  lie  at  the 
same  approximate  level  as  the  mooring  values.  In  Figure  8, 

we  show7  scatterplots  of  the  two  measures  of  model  |u/,|3 
values  at  the  moonng  sites  plotted  against  the  mooring 
values.  As  in  Figure  4,  the  scatter  is  considerable,  the 
number  of  low  outliers  exceeds  the  number  of  high  outliers, 
and  tighter  fits  oeeur  when  the  best  model  values  are  used, 
especially  in  the  ease  of  global  POP. 

[34]  Onee  again,  we  eompute  values  of  7,  using  both 
interpolated  and  best  values,  but  this  time  aeross  the  382 
global  moonng  sites,  and  using  only  the  three  global 
simulations  (Table  4).  In  shallow  and  intermediate  depths, 
7  is  small  for  all  three  global  models,  except  when  global 
POP  is  measured  by  its  best  values.  In  abyssal  waters,  7 


Table  3.  Values  of  7  Computed  in  Subregions  of  the  North 
Atlantic3 


Region 

DANLOM 

NANLOM 

Global 

POP 

North 

Atlantic 

POP 

East  Greenland  Labrador  Sea 

0.051 

0,057 

0.89 

0.99 

Eastern  North  Atlantic 

0.19 

0.058 

0.99 

0.91 

Gulf  Stream 

0.84 

0.68 

0.92 

11 

South  America 

0.31 

0.31 

0.30 

0.39 

a Values  were  computed  from  best  |u/>|3  values  of  the  simulations  in  the 
four  subregions  of  the  North  Atlantic  domain  shown  in  Figure  6. 
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lb)  Best  model  values 


Seafloor  depth  at  mooring  sites  (m) 


Figure  7.  As  in  Figure  3  except  that  the  comparison  is 
global.  Only  the  three  global  simulations  are  shown.  All  382 
near-bottom  current  meters  in  the  global  database  are  used. 


takes  on  order  one  values  (from  0.37  to  0.73)  for  all  three 

simulations  and  for  both  measures  of  model  |u/>|3  values. 
The  best  random  values  of  global  POP  (sec  the  last  line  in 
Table  4)  arc  considerably  lower  than  the  best  values,  for  all 
three  depth  ranges.  The  abyssal  best  random  values  for 
DANLOM  and  NANLOM  are  0.46  and  0.25,  respectively, 
lower  than  the  interpolated  and  best  values  in  the  abyss  for 
these  models.  The  DANLOM  North  Atlantic  abyssal  best 
random  7  value  is  the  only  one  of  several  best  random 
values  we  have  calculated  that  is  not  smaller  than  the 
corresponding  best  and  interpolated  7  values.  All  of  this 
suggests  to  us  that  the  simulations  used  in  this  paper  do 

have  some  skill  in  putting  large  |u/>|3  values  in  approxi¬ 
mately  correct  locations.  To  add  to  this  point,  we  reiterate 
the  results  of  Hurl  hurt  and  Hogan  [2000]  and  Hogan  and 
Hurlburt  [2000],  that  abyssal  eddy  kinetic  energies  are  a 
strong  function  of  horizontal  resolution.  A  glance  at 
Figure  12  of  Hurlburt  and  Hogan  [2000]  makes  clear  that 
if  7  values  were  to  be  calculated  from  NLOM  simulations 
run  with  horizontal  resolutions  of  1/8°  or  less,  they  would 
be  much  lower  than  the  7  values  calculated  here.  This  also 
suggests  that  the  simulations  we  have  chosen  to  analyze  in 


this  paper  have  some  degree  of  skill,  more  so  than  other 
simulations  we  could  have  chosen. 

[35]  The  globally  integrated  dissipations  [D]  for  the  three 
global  simulations  are  shown  in  Table  5.  As  in  the  North 
Atlantic,  most  of  the  dissipation  in  DANLOM  and  NAN¬ 
LOM  takes  place  in  abyssal  waters,  while  about  70%  of  the 
dissipation  in  global  POP  takes  place  in  shallow  and 
intermediate  waters.  Depend  ingon  whether  native  or  com¬ 
mon  values  of  cd  are  used,  the  [D]  values  integrated  over  all 
depths  range  from  0.14-0.65  TW  in  the  three  simulations. 

[36]  Globally  integrated  dissipations  broken  down  by 
latitude  instead  of  by  seafloor  depth  arc  displayed  in 
Table  6  In  NANLOM,  the  dissipation  in  tropical  latitudes 


Figure  8.  As  in  Figure  4  except  that  the  comparison  is 
global.  Only  the  three  global  simulations  arc  shown  All  382 
near-bottom  current  meters  in  the  global  database  are  used 
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Table  4.  As  in  Table  1,  but  Computed  Across  All  382  Near- 
Bottom  Current  Meters  in  the  Global  Database3 


Model 

Measure 

1 

Shallow 

Intermediate 

Abyssal 

DANLOM 

Interpolated 

0.0041 

0.031 

0.73 

DANLOM 

Best 

0.0037 

0.060 

0.62 

NANLOM 

Interpolated 

0.0043 

0.028 

0.37 

NANLOM 

Best 

0.0023 

0.060 

0.40 

Global  POP 

Interpolated 

0.045 

0.15 

0.40 

Global  POP 

Best 

0.52 

0.83 

0.46 

Global  POP 

Best  random 

0.048 

0.10 

0.15 

aThe  number  of  moorings  M  equals  14,  94,  and  274  for  shallow, 
intermediate,  and  abyssal  depths,  respectively. 


is  about  14%  of  the  total.  In  global  POP  and  DANLOM,  the 
tropical  dissipation  represents  about  30%  and  40%,  respec¬ 
tively,  of  the  total.  The  large  tropical  dissipations  occur  in 
different  locations  in  global  POP  than  in  DANLOM.  As 
discussed  previously,  DANLOM,  unlike  the  other  simula¬ 
tions,  has  large  patches  of  high  dissipation  in  abyssal  near- 
equatorial  regions.  In  global  POP,  the  tropical  dissipations 
arc  largest  in  regions  of  strong  currents  along  continental 
edges  (sec  Figure  9,  to  be  described  shortly).  In  DANLOM 
and  global  POP,  about  20%  of  the  dissipation  occurs  in  the 
Southern  Ocean,  while  in  NANLOM,  about  40%  of  the 
dissipation  occurs  there.  Similarly  substantial  ratios  of 
Southern  Ocean  to  total  dissipation  were  seen  in  the 
observational ly  based  estimates  of  Sen  el  al.  [2008].  How¬ 
ever,  as  in  the  paper  by  Sen  et  al.  [2008],  none  of  the 
Southern  Ocean  [D]  values  is  as  large  as  the  wind-power 
input  into  the  Southern  Ocean,  which  represents  about  70% 
of  the  global  integral  [Wunsch,  1998]. 

[37]  Global  maps  of  the  dissipation  rates  computed  from 
the  three  global  simulations,  and  from  the  [D^]  observa- 
tionally  based  estimate  of  Sen  et  al.  [2008]  (adapted  here 
from  their  Figure  3a),  arc  shown  in  Figure  9.  A  drag 
coefficient  value  of  q_=  0  0025  is  used  to  construct  all 
four  maps.  Note  that  [£>3]  =  0.21  TW.  Some  of  the  well- 
known  abyssal  currents,  such  as  the  Gulf  Stream,  Kuroshio, 
Agulhas,  and  Antarctic  Circumpolar  Current,  arc  clearly 
visible  in  all  four  maps.  The  deep  currents  arc  generally 
more  pronounced  (in  both  areal  extent  and  intensity)  in  the 
DANLOM  map  than  in  the  other  three  maps.  On  the  other 
hand,  Figure  9c  shows  several  places  where  global  POP  has 
strong  currents  in  shallow  waters,  which  are  not  present  in 
DANLOM  or  NANLOM.  These  arc  not  present  in  the 
obscrvationally  based  estimate  cither,  since  this  only  covers 
locations  with  seafloor  depths  of  3000  m  or  greater.  The 
DANLOM  map  shows  many  patches  of  high  dissipation  in 
basin  interiors  (sec,  for  instance,  the  patch  located  to  the 
northeast  of  New  Zealand)  and  just  off  the  equator  (espe¬ 
cially  in  the  Atlantic)  which  arc  not  seen  in  the  other  maps. 
These  features  arc  in  part  responsible  for  the  larger  global 
dissipation  rate  and  best  random  7  values  computed  from 
DANLOM.  The  visual  agreement  between  model-based  and 
obscrvationally  based  estimates  (or  amongst  different  model 
estimates)  does  not  appear  to  be  as  close  for  abyssal  flows 
as  it  is  for  surface  flows  (compare  Figure  9  here  to  Figure  12 
from  Maltrud  and  McClean  [2005]).  Nevertheless  there  are 
enough  similarities  between  the  maps  shown  in  Figure  9  to 
give  us  some  degree  of  confidence  in  the  result  that  the 
globally  integrated  dissipation  rates  arc  substantial. 


[38]  The  area  of  low  dissipation  in  the  southeast  Pacific  of 
DANLOM  results  from  an  intentional  suppression  of  the 
abyssal  velocities  in  order  to  prevent  an  unrealistic  Antarctic 
Circumpolar  Current  pathway  in  this  area  (sec  Figure  9b). 
This  suppression  was  done  in  the  model  run,  not  in  the 
analysis,  and  thus  cannot  be  undone.  We  estimate  the  impact 
of  this  patch  of  suppressed  velocities  on  the  DANLOM  [ D ] 
values  in  two  ways:  by  computing  the  ratio  of  the  area  of  the 
patch  to  the  area  of  the  remaining  grid  points  in  the  model, 
and  by  computing  the  [D]  value  in  NANLOM  with  veloc¬ 
ities  in  this  patch  set  to  zero.  Both  methods  suggest  that  the 
patch  reduces  the  DANLOM  [D]  value  by  about  3%. 

[39]  As  in  section  5.1,  we  conclude  with  an  analysis  of  a 
few  subregions  having  strong  currents.  Figure  1 0  shows  the 

locations  of  current  meters  for  which  0.0025  |u/,  | 1  exceeds 

1  mW  m~  2  in  regions  we  denote  by  ‘‘Kuroshio”,  “Southern 
South  Atlantic”,  and  “Southern  Africa”.  Table  7  shows  the 
best  7  values  in  these  regions  computed  from  the  three 
global  simulations.  The  7  values  range  from  0.12  to  0.74, 
with  DANLOM  7  values  being  near  order  one  most 
consistently.  For  global  POP,  interpolated  7  values  (not 
shown)  arc  considerably  lower  (by  factors  of  5-7)  than 
the  best  7  values.  Consistent  with  appearances  in  Figure  9, 
the  7  values  suggest  that  some  prominent  abyssal  currents 
arc  too  weak  in  NANLOM  and  in  global  POP 

6.  Frequency  Content  of  Model  Flows  and 
Effects  of  Infrequent  Temporal  Sampling 

[40]  The  calculations  presented  thus  far  arc  based  on 
model  velocities  that  are  not  low-pass  filtered  The  large 
size  of  the  model  outputs  makes  computing  low-pass 
filtered  velocities  at  all  of  the  model  grid  points  impractical. 
However,  we  did  test  the  frequency  content  of  the  velocities 
at  a  small  number  of  locations,  by  writing  out  NANLOM 
velocities  at  the  abyssal  mooring  sites,  and  applying  a  15- 
day  low-pass  filter  before  computing  the  |u^|3  values.  (We 
cannot  apply  a  3-day  low-pass  filter,  as  we  did  with  the 
current  meter  data,  since  the  model  output  is  every  3  days). 

The  average  of  |u^|3  taken  over  all  the  deep  water  mooring 
sites  is  reduced  by  only  6%  from  the  values  presented 
earlier,  demonstrating  that  the  abyssal  NANLOM  bottom 
layer  flows,  at  least,  arc  predominantly  lowr  frequency.  In 
contrast,  if  we  subsamplc  current  meter  data  every  3  days 
(the  NANLOM  sampling  period)  and  then  apply  a  15-day 

filter,  the  average  of  |u/,|3  computed  over  all  of  the  mooring 
sites  is  reduced  by  much  larger  amounts,  up  to  40%, 
depending  on  the  depth  ranges  considered  Alternatively, 
if  we  use  the  mooring  data  as  is  (sampled  very  frequently) 


Table  5.  As  in  Table  2,  but  Computed  Globally  for  the  Three 
Global  Simulations 


Model 

Cd 

[D] 

Shallow 

Intermediate 

Abyssal 

Total 

DANLOM 

Native 

0.0025 

0.041 

0.60 

0.65 

NANLOM 

Native 

0.0017 

0.017 

0.14 

0.16 

Global  POP 

Native 

0.085 

0.018 

0.042 

0.14 

DANLOM 

Common 

00021 

0.034 

0.50 

0.54 

NANLOM 

Common 

0.0022 

0.021 

0.18 

0.20 

Global  POP 

Common 

0.17 

0.036 

0.085 

0.29 
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Figure  9.  Maps  of  dissipation  rateJmW  m-2)  in  the  three  global  models  (a)  DANLOM,  (b)  NANLOM, 
and  (c)  global  POP,  and  (d)  in  the  [D^]  observational  ly  based  estimate  of  Sen  et  al  [2008].  Figures  9a,  9b, 
9e,  and  9d  extend  to  65°N.  A  common  value  of  cd-  0.0025  was  used  to  make  Figures  9a,  9b,  9e,  and  9d. 
Note  that  the  Sen  et  al  [2008]  estimate  shown  here  only  covers  seafloor  depths  exceeding  3000  m. 
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Table  6.  As  in  Table  5,  but  Divisions  Are  by  Latitude  Instead  of 
by  Seafloor  Depth3 


m 


Model 

cd 

Tropical 

Extratropical 

Southern  Ocean 

Total 

DANLOM 

Native 

0.27 

0.38 

0.13 

0.65 

NANLOM 

Native 

0.022 

0.14 

0.064 

0.16 

Global  POP 

Native 

0.043 

0.10 

0.030 

0.14 

DANLOM 

Common 

0.22 

0.32 

0.11 

0.54 

NANLOM 

Common 

0.028 

0.17 

0.080 

0.20 

Global  POP 

Common 

0.087 

0.21 

0.062 

0.29 

“Tropical  denotes  latitudes  equatorward  of  23.5°,  extratropical  denotes 
latitudes  poleward  of  23. 5°,  and  Southern  Ocean  denotes  latitudes  south  of 
40°S.  Total  denotes  integration  over  all  latitudes. 


and  low-pass  it  with  just  a  3 -day  filter,  the  average  of  |u/,|3 
calculated  over  all  of  the  mooring  sites  is  reduced  by  20% 
or  more  [Sen  et  al,  2008].  Both  exercises  demonstrate  that 
high-frequency  motions  arc  much  more  prominent  in  cur¬ 
rent  meter  data  than  in  the  simulations  used  in  this  paper. 

[41]  The  interannual  variability  of  globally  integrated 
dissipation  appears  to  be  small.  We  analyzed  year  2002  of 
global  POP,  and  found  little  difference  in  the  dissipation 
rates  from  those  presented  in  section  5  using  year  2003. 
Analysis  of  a  separate  year  of  DANLOM  output  from  that 
used  here  indicates  that  it  also  shows  little  interannual 
variability  in  the  dissipation  rate. 


(a)  Kuroshio  (b)  Southern  South  Atlantic 


(c)  Southern  Afnca 
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Figure  10.  As  in  Figure  6,  but  for  three  subregions  outside 
of  the  North  Atlantic  domain,  (a)  Kuroshio  (Af  =  8,  Daverage  = 
1 2).  (b)  Southern  South  Atlantic  (M  -  22,  Daverage  =  22). 
(c)  Southern  Africa  (M  =  13,  Daverage  =  3.5). 


Table  7.  As  in  Table  3,  but  for  the  Three  Subregions  Outside  of 


the  North  Atlantic  Shown  in  Figure  10 


Region 

DANLOM 

NANLOM 

Global  POP 

Kuroshio 

0.40 

0.19 

0.32 

Southern  South  Atlantic 

0.31 

0.16 

0.12 

Southern  Africa 

0.71 

0.74 

0.36 

[42]  High-frequency  motions  will  be  aliased  into  the 
model  output,  because  of  the  infrequent  temporal  sampling 
of  model  snapshots.  However,  McClean  et  al.  [2002] 
showed  that  the  error  introduced  by  infrequent  sampling 
was  minimal,  for  the  ease  of  surface  kinetic  energy  in  daily 
snapshots  of  North  Atlantic  POP.  Here  we  examine  the 
impact  of  subsampling  on  one  year  of  output  from  another 
simulation  of  nonassimilative  NLOM,  which  we  denote  by 
NANLOM2.  The  NANLOM2  run  is  similar  to  the  NAN- 
LOM  run,  except  that  it  contains  the  same  patch  of  sup¬ 
pressed  abyssal  velocities  in  the  South  Pacific  seen  in 
DANLOM  (whieh  is  why  we  didn’t  use  it  in  section  5), 
has  a  slightly  different  Ca  value  (0.003),  and  most  impor¬ 
tantly,  saved  snapshots  of  model  output  more  frequently, 
every  6  hours.  For  the  sake  of  comparison,  we  also  examine 
the  impact  of  subsampling  on  the  unfiltered  current  meter 
data,  whieh  as  we  have  discussed  contains  stronger  high- 
frequency  flows  than  the  models. 

[43]  Figure  1 1  shows  the  main  results  of  our  investiga¬ 
tions,  histograms  of  the  quantity 


m;2  -  iu„i;, 

W/, 


(4) 


where  72  represents  a  coarser  sampling  interval,  and  71  a 
finer  interval,  of  the  same  data  set.  Figure  11a  shows 
histograms  of  a  for  the  data  set  of  near-bottom  eurrent 
meters,  where  72  =  6  hours  and  71  represents  the  ” source” 
sampling.  The  souree  sampling  is  the  sampling  given  in  the 
mooring  database,  which  ranges  from  15  min  to  3  hours 
depending  on  the  particular  mooring.  To  avoid  small 
numbers  in  the  denominator  of  a,  we  also  use  only  those 
moorings  for  which  0.0025  p  exeeeds  1  mW  m 
There  arc  224  such  moorings.  Six  of  the  moorings  arc 
sampled  at  3-hour  intervals,  while  the  remainder  arc 
sampled  at  intervals  of  2  hours  or  finer,  most  commonly 
1  hour.  Figure  1 1  a  shows  that  for  most  of  the  current  meters, 
subsampling  at  6-hour  intervals  yields  nearly  the  same  value 

of  | u?, | 3  as  obtained  from  the  more  frequent  source 
sampling.  This  suggests  that  the  6-hour  sampling  interval 
in  NANLOM2  is  sufficiently  fine  for  it  to  serve  as  an 
adequate  “control”  (assumed  to  have  relatively  little 
aliasing)  for  an  analysis  on  the  impacts  of  infrequent 
sampling.  The  histogram  of  a  values  computed  from  daily 
subsamplcs  (72  —  1  day)  and  71  still  representing  the  source 
sampling  is  notieably  wider  (Figure  1 1  b).  Therefore  a 
substantial  number  of  moorings  are  adversely  affected  by 
1-day  subsampling.  Because  the  difference  between  6  hourly 
subsampling  and  souree  sampling  is  small  (Figure  1  la),  and 
because  NANLOM2  uses  6  hourly  sampling,  we  use  the 
6  hourly  sampling  as  a  baseline  (71  =  6  hours)  in  the  mooring 
data  to  test  the  effects  of  the  infrequent  temporal  sampling  of 
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Figure  11.  Histograms  of  a  =  where  72  represents  a  coarser  sampling  interval  and  71 

fud/i 

represents  a  finer  interval  of  the  same  data  set  Results  shown  (a-c)  for  moored  current  meters  and  (f-j)  for 
model  grid  points  in  NANLOM2  (a  6  hourly  sampled  nonassimilativc  simulation  of  NLOM).  In 
Figures  11a,  lib,  lie,  lid,  and  lie,  source  denotes  the  original  mooring  data,  sampled  at  intervals  from 
1 5  min  to  3  hours  depending  on  the  mooring.  All  other  sampling  intervals  are  given  at  the  top  of  Figures  1 1  a, 
1  lb,  1  le,  1  Id,  1  le,  Ilf,  1  lg,  1  lh,  Hi,  and  1  lj.  Only  the  224  moorings  and  the  2094302  model  grid  points  for 
which  0.0025  p  |ua|;|  exceeds  1  mW  m-2  arc  used. 


DANLOM,  NANLOM,  and  POP.  In  Figures  lie,  lid,  and 

I  le  we  set  72  =  1,3,  and  10  days,  the  sampling  periods  of 
these  three  simulations,  respectively.  The  daily  minus  6 
hourly  histogram  (Figure  lie)  looks  very  similar  to  the  daily 
minus  source  histogram  (Figure  lib),  as  expected.  The 
histograms  made  with  72  =  3  and  10  days  (Figures  1  Id  and 

I I  c,  respectively)  become  progressively  wider,  meaning  that 
the  \Ub[  values  at  individual  current  meters  become  less  and 


less  reliable  as  the  sampling  period  increases.  Note  that  the 
results  plotted  in  Figures  11a,  lib,  lie,  lid,  and  lie  arc 
dependent  on  the  record  lengths,  which  vary  between  current 
meters,  as  well  as  on  the  subsampling  period. 

[44]  In  Figures  Ilf,  1  lg,  1  lh,  Hi,  and  1  lj,  a  is  calculated 
from  NANLOM2,  using  71  =  6  hours  and  72  =  1,  3,  10,  20, 
and  40  days,  respectively.  As  in  the  mooring  analyses,  we 
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only  use  model  grid  points  for  which  0.0025  p  |u^|yj 
exceeds  1  mW  m  2.  There  are  2094302  such  model  grid 
points.  For  12  =  1  and  3  days,  the  histograms  of  a 
(Figures  Ilf  and  llg,  respectively)  are  quite  narrow,  dem¬ 
onstrating  that  the  subsamplings  used  in  DANLOM  and 
NANLOM  are  adequate  for  our  purposes.  For  12  =  10  days, 
the  sampling  period  for  POP,  the  histogram  (Figure  1  lh)  has 
roughly  the  same  width  as  the  histogram  of  the  daily 
subsamplcd  mooring  data  (Figures  lib  and  lie).  The 
histograms  for  20-  and  40-day  subsampling  of  NANLOM2 
(Figures  1 1  i  and  1  lj,  respectively)  are  of  comparable  widths 
to  the  histograms  for  3-  and  10-day  sampling,  respectively, 
of  the  current  meter  data.  All  of  this  further  suggests  that 
these  simulations  do  not  have  as  much  high-frequency 
variability  as  the  current  meters,  making  the  aliasing  prob¬ 
lem  less  severe  than  it  is  in  the  current  meters  for  the  same 
sampling  interval. 

[45]  The  histograms  of  both  current  meter  data  and  model 
output  display  approximately  similar  numbers  of  negative 
and  positive  a  values,  meaning  that  when  subsampling 

matters  it  can  either  decrease  or  increase  the  |u^|3  estimates. 
The  average  over  all  of  the  moorings  (or  model  grid  points) 
of  | Ub | 3  is  therefore  relatively  unaffected.  In  the  analyses 
shown  in  Figures  11a,  lib,  lie,  lid,  and  lie,  the  averages 
over  the  subsampled  (12)  mooring  data  are  within  8%  or  less 
of  the  averages  computed  from  the  better  sampled  (71)  data. 
If  we  use  only  abyssal  moorings  the  comparison  is  even 
closer.  In  the  analyses  shown  in  Figures  Ilf,  llg,  llh,  lli, 
and  1  lj,  the  area-weighted  averages  over  all  of  the  NAN- 

LOM2  grid  points  of  the  subsampled  |u/>|3  values  lie  within 
3%  or  less  of  the  area-weighted  averages  of  the  6  hourly 
values.  The  latter  exercise  demonstrates  that  the  [D]  esti¬ 
mates  given  in  section  5  are  not  severely  affected  by  the 
infrequent  temporal  sampling  of  the  models. 

7.  Wind-Power  Input 

[46]  To  test  whether  the  energy  inputs  into  the  models  are 
of  the  expected  magnitude,  we  calculate  the  wind-power 
input  for  the  NANLOM  simulation.  We  use  the  formula 

W  —  J  r  u gdA,  (5) 

where  r  is  wind  stress  and  is  the  gcostrophic  velocity  at 
the  sea  surface,  calculated  from  centcrcd-diffcrcncing  of  the 
model  sea  surface  height  fields.  The  latitude  band  3°S-3°N 
is  excluded  from  the  calculation,  since  the  gcostrophic 
relation  is  not  valid  on  the  equator.  The  average  wind  work 
over  year  2002  calculated  from  NANLOM  snapshots  put 
out  every  3  days  is  W  =  0  61  TW,  comparable  to  but 
somewhat  lower  than  the  0.75-0.9  TW  estimated  by  Scott 
and  Xu  [2008]  for  the  wind  work  poleward  of  3°. 

8.  Summary  and  Discussion 

[47]  Computations  on  the  output  of  three  existing  global 
high-resolution  simulations  strongly  suggest  that  the  energy 
dissipation  rate  of  the  oceanic  general  circulation  by  qua¬ 
dratic  bottom  boundary  layer  drag  represents  a  significant 
fraction  of  the  ~1  TW  wind-power  input  transformed  into 


gcostrophic  flows.  The  range  of  dissipation  estimates  over 
all  water  depths  obtained  from  the  simulations  is  0.14- 
0.65  TW,  quite  comparable  to  the  0.21-0.83  TW  of 
dissipation  in  waters  deeper  than  3000  m  estimated  by 
Sen  et  al  [2008]  from  observations  alone.  In  addition,  the 
patterns  of  dissipation  estimates  obtained  from  the  models 
arc  qualitatively  similar  to  those  obtained  from  observations 
in  the  paper  by  Sen  et  al  [2008].  Both  the  model-based  and 
observational ly  based  maps  show  dissipation  concentrated 
in  strong  flows  such  as  the  Gulf  Stream,  fCuroshio,  Antarctic 
Circumpolar  Current,  etc. 

[48]  Prior  to  the  computations  of  dissipation,  we  com¬ 
pared  the  time  average  of  the  cube  of  bottom  velocities 

(|u/,|3)  computed  from  the  simulations,  with  |u/,|3  values 
computed  from  near-bottom  moored  current  meters.  The 

comparisons  focus  on  \u^\  bceause  it  is  proportional  to  the 
dissipation  rate  The  comparisons  indicate  that  numerical 
simulations  generally  have  difficulty  matching  individual 
current  meter  records  However,  the  ratio  of  the  average  of 

the  model  |u/,|3  values  at  274  moored  near-bottom  current 
meter  sites  in  the  abyssal  ocean  to  the  average  computed 
from  the  current  meters  themselves  is  of  order  one  (0.37- 
0.73)  for  all  three  simulations.  This  suggests  that  the  model 
estimates  of  areally  integrated  dissipation  rate  should  be 
correct  to  within  an  order  of  magnitude. 

[49]  Since  the  ratio  just  described  is  less  than  one  in  the 
abyss,  and  is  much  less  than  one  in  shallow  waters  for  the 
NLOM  simulations,  some  of  our  computed  dissipation  rates 
are  likely  to  be  underestimates.  The  wind-power  input  into 
NANLOM  was  found  to  be  M).  15-0.3  TW  lower  than 
observational  estimates  [Scott  and  Xu ,  2008],  further  rein¬ 
forcing  the  possibility  that  some  of  our  computed  dissipa¬ 
tion  rates  are  likely  to  be  underestimates.  On  the  other  hand, 
there  arc  two  reasons  to  believe  that  DANLOM  in  particular 
may  overestimate  the  dissipation  DANLOM  shows  dissi¬ 
pation  in  regions  where  it  is  not  found  in  the  other  models. 
In  addition,  the  data  assimilation  methodology  used  in 
DANLOM  is  not  energy  conserving,  such  that  the  data 
assimilation  process  constitutes  an  extra,  unphysical  energy 
source  for  DANLOM. 

[50]  The  estimated  dissipation  rates  in  the  Southern 
Ocean  represent  a  substantial  fraction  of  the  total  dissipa¬ 
tion.  However,  they  fall  far  short  of  the  wind-power  input 
into  that  region,  suggesting  that  cither  bottom  drag  in  other 
regions  of  the  ocean,  or  mechanisms  other  than  bottom  drag 
acting  in  the  Southern  Ocean,  balance  a  substantial  fraction 
of  the  local  input.  Similarly,  since  the  lower  end  of  our 
global  dissipation  estimates  balances  only  about  15-20%  of 
the  global  wind-power  input,  we  conclude  as  we  did  in  the 
paper  by  Sen  et  al  [2008]  that  dissipation  mechanisms  other 
than  bottom  boundary  layer  drag  arc  also  likely  to  be 
significant  contributors  to  the  global  energy  budget  of  the 
general  circulation.  In  the  models  utilized  here,  energy 
which  is  not  dissipated  by  bottom  drag  is  dissipated  at 
small  horizontal  scales,  cither  by  eddy  viscosity  (NLOM), 
or  biharmonic  viscosity  (POP).  The  extent  to  which  the 
transfer  of  gcostrophic  energy  to  smaller  scales  represents  a 
significant  dissipator  in  the  real  ocean  is  an  area  of  active 
research  [c.g.,  Muller  et  al ,  2005;  Polzin ,  2008].  It  is  quite 
possible  that  small-scale  damping  operators  such  as  eddy 
viscosity  will  ultimately  be  seen  as  extremely  crude  proxies 
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for  the  actual  physics  of  such  transfers.  The  generation 
of  internal  waves  by  geostrophie  flows  over  rough  topog* 
raphy  [e.g.,  Naveira-Garabato  et  ai ,  2004;  Marshall  and 
Naveira-Garabato ,  2008;  Nikurashin ,  2008]  is  another 
potentially  important  dissipation  mechanism  that  is  not 
included  even  in  a  crude  way  in  the  models  utilized  in  this 
study. 

[5 1 ]  High-resolution  ocean  models  continue  to  improve 
and  as  they  do  so  the  computations  performed  here  will  bear 
repeating.  As  discussed  in  this  paper,  higher  resolution  in 
both  the  vertical  and  horizontal  directions  will  almost 
certainly  improve  the  comparison  of  bottom  flows  in 
models  to  those  in  observations.  Improvements  in  bottom 
bathymetry  products  and  in  data  assimilation  schemes,  and 
inclusion  of  effective  parameterizations  of  other  dissipation 
mechanisms  in  addition  to  bottom  boundary  layer  drag,  are 
also  likely  to  have  a  positive  impact  in  this  regard. 
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